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so: 

r—{ t We calculate na phase-shifts and scattering observables in the resonating-group method, 

using the nuclear-matter G-matrix of an SUe quark-model NN interaction. The G-matrix 
is generated in the recent energy-independent procedure of the quark-model AW interaction 
' with the continuous prescription for intermediate spectra, by assuming an appropriate Fermi 

momentum kF = 1.2 fm -1 . The na RGM interaction kernels are evaluated with explicit 
treatments of the nonlocality and momentum dependence of partial-wave G-matrix compo- 
i^h ■ nents. The momentum dependence of the G-matrix components is different for each of the 

I ' nucleon-exchange and interaction types. Without introducing any artificial parameters ex- 

cept for kF, the central and spin-orbit components of the na Born kernel are found to have 
reasonable strengths under the assumption of a rigid translationally invariant shell-model 
wave function of the a-cluster. The characteristic behaviors of three different exchange terms, 
1 1 1 corresponding to knockout, heavy-particle pickup and nucleon-rearrangement processes, are 

essentially the same between the case of previous local effective NN forces and the case of 
C^A ' nonlocal G-matrix NN interactions. 

> ' 

o : 

, §1. Introduction 

(N ■ 

Microscopic cluster models have been successfully used to describe structure 
and reactions of light nuclear systems. For instance, low-energy na scattering is well 
described by the na resonating-group method (RGM) with various model spaces 
and effective nucleon-nucleon (NN) forces. These effective forces usually incorpo- 
rate only the central and spin-orbit (LS) forces. The very strong one-pion exchange 
tensor force of the bare NN interaction is renormalized to the S E central force ef- 
fectively. The usage of these effective forces is justified largely by the success of 
the RGM calculations, in which the existence of ample experimental data for the 
NN scattering and light nuclear systems are prerequisite. On the other hand, the 
experimental data in the hypernuclear systems are still not yet sufficient and basic 
baryon-baryon interactions are not well known because of the technical difficulties of 
strangeness experiments. From the theoretical side, some progress is made with the 
study of baryon-baryon interactions and with accurate calculational techniques to 
solve few-body systems. One of these attempts is our effort to construct a unified set 
of quark-model Bg-Bg potentials 1 ) and to apply them to solve few-body systems. 2 ) 
We now need a procedure to link bare and effective interactions through some ef- 
fective interaction theory such as the G-matrix formalism. In this paper we present 
such a calculational scheme for light nuclei, in which we introduce no intermediate 
localized effective potential, and directly use G-matrices in nuclear matter. First, we 
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reexamine the well-studied problem of na RGM calculations from the viewpoint of 
our method to establish the reliability to proceed to hyperon-nucleus systems. The 
latter applications are to be reported in a successive paper. 

We will carry out na RGM calculation using a quark-model G-matrix NN inter- 
action and a framework that has recently been developed for a-cluster folding in the 
study of the baryon-octet (Bg) a interaction. 3 ) In this framework, the partial- wave 
components of the G-matrix are explicitly used to generate the direct and knock-o 
terms of the na RGM kernel, without making any kind of local approximation for 
the G-matrix. The center-of-mass (cm.) motion of two interacting nucleons in the 
na system is correctly treated for the Galilean non-invariant G-matrix interaction. 
The G-matrix and momentum-dependent single-particle (s.p.) potentials are pre- 
determined by solving the Bethe-Goldstone equation in symmetric nuclear matter. 4 ) 
The only assumption is a constant Fermi momentum hp for generating the G-matrix 
interaction that can be used as an effective interaction in light nuclear systems. We 
use the energy-independent version of the quark-model NN interaction to calculate 
the G-matrix, but the difference from the previous energy-dependent version in Ref. 
4) is very little. The original quark-model NN interactions, fss2 and FSS, used in 
this paper are a unified model for full octet-baryons, 1 ) which have achieved accurate 
descriptions of the NN and YN scattering observables. In particular, the NN in- 
teraction of the most recent model fss2 5 ) is sufficiently accurate in comparison with 
those of modern realistic meson-exchange models. 

Since the na RGM kernel involves various nucleon-exchange terms, we need to 
extend the previous folding formula, starting from the transformation formula de- 
veloped in Appendix A of Ref. 6). In general, the interaction kernel involves five 
different interaction types for two-cluster configurations with a common width pa- 
rameter of (Os) clusters. Among them, the exchange terms, called the 15" and IS" 
types in this paper, do not appear in the previous hyperon-a interaction and need 
a special treatment in the present na study. These terms correspond to the so- 
called heavy-particle pickup process 7 -* and play an important role in the backward 
increase in the differential cross sections in the low-energy and intermediate-energy 
regions. On the other hand, the knock-on term specified by the ID- type gives 
an energy-dependent extra attraction to the direct potential specified by 0D + , and 
contributes mainly to the forward direction in the usual RGM treatment with effec- 
tive ./ViV forces. In the present formalism, all these interaction kernels have their 
own momentum dependence for the Galilean non-invariant two-nucleon interaction. 
We can calculate the momentum dependence explicitly in the analytic form. A 
very important starting-energy dependence of the G-matrix is renormalized into this 
momentum dependence and the dependence to the relative momentum, and to the 
Fermi momentum hp as well. Owing to the explicit evaluation of all the interaction 
kernels, the existence of the Pauli forbidden (Os) state between the neutron and 
the a cluster is strictly preserved. We will find that the present procedure gives 
rather reasonable strengths of the na interaction for both the central and LS com- 



*' We will use in this paper the terminology "knock-on term" to specify the interaction terms 
responsible for the "knockout process" in the nuclear reaction mechanism discussed in Ref. 7) . 
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ponents, and reproduces reasonably well the empirical S'-wave and P-wave phase 
shifts and the low-energy differential cross sections and polarizations below the neu- 
tron incident energy E n ~ 30 MeV. For higher energies, we investigate the na Born 
kernel and find that the characteristic behavior of three different groups of exchange 
terms, corresponding to knockout, heavy-particle pickup and nucleon-rearrangement 
processes, found in previous studies using local effective NN forces, 7 ) is essentially 
unchanged even in the present study with the nonlocal G-matrix interaction. 

The na RGM has been examined by many authors from many different view- 
points. In fact, early studies have naturally paid full attention to the adequacy of 
the model space and the NN force used in the calculation. For example, Sugie et 
a/. 8 ) studied the contribution of tensor force, trying to explain the large energy 
splitting of the J w = 3/2 _ and l/2~ states. Since the inherent LS force of the NN 
interaction has been found after this calculation, the authors did not include it and 
found that tensor force can account for only about 30% of the observed splitting. 
Later calculations by Kanada et al. 9 ^ and Omojola 10 ) included the LS force and 
have found that it is important to deal with the D-state components of the a cluster 
and realistic NN interactions. In these calculations, an approximate Hamada and 
Johnston potential is used in the Gaussian form. They reproduced the low-energy 
na phase shifts reasonably well. The channel coupling effect of the na and d 3 H 
configurations are examined by Heiss and Hackenbloich. 11 ) Thompson and Tang 7 ) 
analyzed the properties of the na RGM exchange kernels for an effective central NN 
force, and classified them into three different groups of terms, called the knockout, 
heavy-particle pickup and nucleon-rearrangement terms. Chwieroth et al. l2 > and 
later calculations 13 ^' 14 ) have clarified that the distortion effect of the a cluster also 
has an appreciable effect even in single-channel na RGM calculation. Nevertheless, 
the effect is rather minor, in comparison with the other two-cluster systems, since 
the a particle is tightly bound. After all of these investigations, it is gradually rec- 
ognized that single-channel na RGM calculation, using a rigid (0s) 4 a-cluster wave 
function and a simple effective NN force of the central and LS types, reproduces the 
low-energy na scattering fairly well, except for some specific energy regions where 
other reaction channels open. This understanding cannot directly be reconciled with 
the recent ab initio calculation by Nollett et al., 15 ^ who claim that the correct P- 
wave spin-orbit splitting of the na scattering phase shifts in the low-energy region 
can only be achieved with the effect of an appropriate three-nucleon interaction. 

The organization of this paper is as follows. In the next section, we first recapit- 
ulate in Sec. 2.1 the standard RGM formulation in the momentum representation, 
together with the G-matrix calculation of the quark-model baryon-baryon interaction 
for symmetric nuclear matter. The full expressions of the na exchange interaction 
kernels for the G-matrix NN interaction are given in Sec. 2.2. The partial-wave 
decomposition is presented in Sec. 2.3. In Sec. 2.4, we discuss the selection of the 
starting energy in the G-matrix calculation. Numerical results of the G-matrix cal- 
culation and the na RGM phase shifts are given in Sec. 3.1. The na scattering cross 
sections and polarization for neutron incident energies less than 30 MeV are com- 
pared with experimental results in Sec. 3.2. The characteristic behavior of the na 
Born amplitudes for higher energies is analyzed in Sec. 3.3. Section 4 is devoted to 
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the summary of this paper. The na RGM Born kernels for a Gaussian-type effective 
NN force are given in Appendix A. In Appendix B, we give the exchange interaction 
kernels of a Galilean non-invariant NN interaction for general systems, composed of 
two (Os)-shell clusters. 

§2. Formulation 

2.1. Lippmann-Schwinger RGM and the G-matrix calculation of symmetric nuclear 
matter 

For the correct treatment of the cm. motion, it is most convenient to formulate 
the RGM in the momentum representation, which we call the Lippmann-Schwinger 
RGM (LS-RGM). 6 ) In this approach, we write the RGM equation in the Schrodinger- 
type equation 

[e-T r -V* GM (e)] X = 0, (2-1) 

where T r is the kinetic-energy operator for the relative motion of two clusters, e = 
E — Ei ni with E- mt being the internal energy, and T/ RGM (e) = Vd + G + eK the 
sum of the direct potential Vb, the exchange kernel G, and the energy-dependent 
term eK, inherent in the RGM formalism. We use the notation K for the exchange 
normalization kernel, which is defined through 

N = (0 in W nt > = 1 + X N M N = 1 - K . (2-2) 

In the na RGM, the spin-isospin factor for the exchange normalization kernel is 
Xn = —1, and Mn is the corresponding spatial part. The total wave function 
is expressed as = A'{(p mt x} using the relative wave function \ hi Eq. (|2-1|) . 
Here, (p mt is the internal wave function and A' the antisymmetrization operator 
between two clusters. For the na system, these are given by 4> mt = <p a and 

2 2 

A' = 1 — Ya=i ~ > 1 — 4P(45) using the internal wave function of the a cluster 
4>a = (t>a h £oo with 0° rb being the spatial part of the translationally invariant (Os) 4 
harmonic-oscillator (h.o.) wave function without the cm. motion. The spin-isospin 
wave functions are denoted by £oo and £i i for the a cluster and the fifth nucleon, 

2 2 

respectively. In this particular case, we have one Pauli-forbidden state |it(o s )) satis- 
fying 

K\u (0s) ) = \u {0s) ) and A'{cf> int u {0s) } = . (2-3) 

The exchange kernel G = G K + G V is composed of the exchange kinetic-energy kernel 
G K and the exchange interaction kernel G v for the total Hamiltonian consisting of 

5 5 

H = J2ti-T G + Y,Vij (2-4) 

i=l i<j 

in the na RGM. Equation (j2- lj) is then solved in the form of the Lippmann-Schwinger 
equation for the half-off shell T-matrix with discretized momentum mesh points, 
which is fully spelled out in Ref. 6). 



na Resonating- Group Calculation with a Quark-Model G -Matrix NN int. 5 



The main task in the LS-RGM is therefore to calculate the plane-wave ma- 
trix elements of various exchange kernels, which we call the RGM Born kernels. 
For the interaction kernels, we separate the two-nucleon force into the spatial and 
spin-isospin parts, using Vij = UijWij. The direct potential Vd and the exchange 
interaction kernel G v are obtained from the general expression 

M(q f , q t ) = ( 6(X G ) e^f r ^ | £ Vij A' | 1 ■ e^ r 

i<j 

= Y, X *T M xr {q f , qi ) . (2-5) 
xT 

Here, the spatial integral M x q-{q^q i ) and the spin-isospin factor X x q- are defined 
by 

M x r(Qf,q l ) = (5{Xg)z x e i( lf r <t> orh \ Uij \l-e i( li- r ct> orb ) , 

r 

x xT = c x ( Zx e i w a i o with (*> j) G r • ( 2 - 6 ) 

i<j 

In the na RGM, (j) orh = 0° rb and £ = £oo £i i ■ These kernel components are specified 
by the number of exchanged nucleons, x = 0, 1, with Cq = 0, C\ = —4, and various 
interaction types T. More specifically, x = with zq = 1 corresponds to the direct 
terms including T = E and D + interaction types and x = 1 with zi = Pus) the 
one-nucleon exchange terms. The interaction types, T = E, S, S', D + and -D_, 
correspond to some specific (i,j) pairs of the nucleons. 16 ^ When an effective NN 
force is used, these Born kernels are most easily calculated using a general trans- 
formation formula developed in Appendix A of Ref. 6). For Gaussian- type effective 
NN forces with the central form 

V (C) = u ^ W = VQ e -«r 2 ^ w + Bp ^ _ HPt _ M P a P T ) , (2-7) 
and the LS form 

V (LS) = u LS( r ) W LS = V LS e ~ K r 2 (yy _ Hp ^ (X . 5) ; (2-8) 

the final results of the na Born kernels are given in Appendix A for completeness. 

The quark-model baryon-baryon interaction is formulated in a similar way to the 
na RGM. 1 ) In this case, the internal wave function is a product of two three-quark 
clusters (3q)-(3q). We solve the G-matrix equation 

G NN (p, q; K, co, k F ) = V$% M (p, q) 

+ (2^ / d k V ^ M ^ k ^e { (k,K- k J) GNN(k,q;K,u,k F ) , (2-9) 

using the energy-independent Born kernel V^ M (p, q) for the (3q)-(3q) system. 4 ) 
Here, V$$ M (p, q) is defined by 17 ) 

V^ M ( P , q) = V D (p, q) + G(p, q) + W(p, q) , (2-10) 



6 



Y. Fujiwara, M. Kohno and Y. Suzuki 



W = -±=(T r + V D + G)-±=-(T r + V D + G) . (2-11) 



with 



Details of this energy-independent treatment of the quark-model baryon-baryon in- 
teraction in the G-matrix formalism will be published elsewhere. In Eq. (|2-9p . we 
use the angle-averaged Pauli operator Q(k,K,kp), and the energy denominator 
e(k, K; uj) is given by 

e(k,K;u;) = u-E N (ki)-E N (k2), (2-12) 

with the starting energy u expressed as 

u = E N ( qi ) + E N (q 2 ). (2-13) 

The s.p. momenta q 1 and q 2 (^l an d ^2) are related to the relative momentum q (k) 
and the cm. momentum K through the conventional relationship q = {q l — q^)/^ 
and K = q 1 + q 2 (k = (fei — k 2 )/2 and K = k\ + k 2 ). The nucleon s.p. potential 
Un{qi) included in the s.p. energy 

EN(qi) = ^qi 2 + U N ( qi ) (2-14) 

is determined self-consistently in the standard procedure. In Eq. (|2-14j) . M^r is the 
nucleon mass. The procedure to include the s.p. potential even in the intermediate 
spectra in Eq. (|2-12j) is called a continuous prescription. 

2.2. Interaction kernels for G-matrix NN interaction 

In this subsection, we derive the interaction Born kernel for the G-matrix NN 
interaction. The starting point is the invariant G-matrix 3 ) expressed aJ*^l 

= \{ \NN\n, I G(p,p'; K,u, k F ) - G(p, -p'; K, uj, k F ) P a P T \ [NN] n , ) 

= gl + g I ss {<Ti-<T 2 ) + h I Q in-(<T l + <T 2 ) + --- . (2-15) 

Here n = \p' xp]/(j?'p sin 9), and the invariant functions g$ (central), gl s (spin-spin), 
/iq (LS), etc. are functions of p = \p\, p' = \p'\, and cos^ = (p -p'), as well as 
the G-matrix parameters K, uj and kp. These are expressed by the partial- wave 
components of the NN G-matrix as in Appendix D of Ref. 6). In the following, 
we will focus on the momentum dependence of the G-matrix, and keep only the 
parameter K in G ! NN (p,p'; K,co,kp), since the explicit dependence depends on the 
interaction type. As in Ref. 3), it is convenient to write the isospin dependence of 
the invariant G-matrix as 

G NN (p,p' ] K) = Gjf^p,p' ] K)^-^ + Gj^(p,p' ] K)^—^ . (2-16) 



*' The invariant NN G-matrix in Eq. (|2-15[) is defined without the factor 2, shown in Eq. (2.4) 
of Ref. 3), since the exchange terms are explicitly calculated in the present na RGM formalism. 
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Table I. Spin-isospin factors X^t f° r the invariant species, Q — (go), ss (g ss ) and LS (ho), in 
Eq. (H7]). The factors for xT = IE, IS and IS' are obtained from = Xis = X?s> = 

(— 1/2)Xq E i . For the LS term, X^ 1 = X^ 1 ■ 2S is assumed with the neutron spin operator 
S. 



I 


Q 


v n i 

A 0E 


v n i 


v n i 







1 


1 


i 

2 


1 


ss 


-3 





3 
2 




LS 




1 


-1 







3 


1 


1 

2 





ss 


3 





3 
2 




LS 




1 


1 



We find it convenient to separate the isospin multiplicity factor (21 + 1) and define 
the spin-isospin factors in Eq. (|2-6j) by 

x£ r = (2I + l)X° T I . (2-17) 

The interaction species f2 = 0, ss and L5", correspond to the invariant functions, go, 
g ss and ho, respectively. Table 1 lists for each of these species. 

The spatial integrals are obtained by assuming a general Galilean non-invariant 
interaction 

(p 1 ,p 2 \u\p[,p' 2 )=S(K-K')-^u(k',q' ] \K\) , (2-18) 

where u is g^, gl s , or h$ in. Here the relative momentum p and the total momentum 
K (and also p' etc. with primes) are related to p 1 and p 2 by p = (p x — p 2 )/2 and 
K = p 1 +p 2 , respectively, and a further transformation from p and p' to k' = p — p' 
and q' = (p + p')/2 is applied. The necessary spatial integrals for the na system 
is obtained from more general expressions given in Appendix B for systems of two 
(Os)-shell clusters. 

For the Galilean non-invariant G-matrix interaction, the subtraction of the 
internal-energy part of the a cluster involves a subtle problem. In the total cm. sys- 
tem of the na system, the O-E-type spatial integral, involving the a-cluster internal- 
energy contribution, becomes momentum-dependent. Similarly, the 1^-exchange- 
type spatial integral also involves a momentum dependence besides the exchange 
normalization kernel M^{qj,q i ). These are explicitly given by 

M 0E (q f , Qj = (27T) 3 5(k) E^ E (q) , M 1E (q f , Qi ) = M N (q f , Qi ) EY E (q) , (2-19) 
where k = qj — q i , q = (q^ + Qj)/2, and E^ E (q) with x = 0, 1 are calculated from 

3 

EYe(i) = {J^y jKe-UK-U) 2 £ { K) , (2-20) 



s 
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with 

3 

{-^y i dpdp ' exp {"^( p2+p/2 )} 9{p > p, > K) - (2 - 2i) 

Here g(p,p'; K) = u(k',q';K) and v is the h.o. width parameter for the (Os) 4 a 
cluster. Owing to the different q dependence involved in the OE and IE types, the 
standard subtraction of the a-cluster internal energy in l/ RGM (e) is not complete 
and is given by 

V RGM (e) = V B + G K + G y + eK , 
Vd = Xqd + Mqo + , 

G v = ^X 0E (Mie - Mis ~ M 1S ^ + X lD _M lD _ , 

Mxe = 2M 0E - Mie - 2M ( ° } (1-K) . (2-22) 

Here we have omitted fl and the isospin sum for simplicity. In deriving Eq. (|2-22l) . 
we have assumed that only the q = part of EQ E (q) contributes to the a internal 
energy and defined the relative energy e using this zero-momentum a energy. In 
Eq. (JEED, M ( § implies 

M^% f , qi ) = (2irf5(k)EV E (0) . (2-23) 

Since the delta function part in M\ E of Eq. (|2-22p is inconvenient for the exchange 
interaction kernel, we assume 

M lE {q f , qi ) = M N (q f , qi ) [2E^ E {q) - E\ E (q)] , (2-24) 

without violating the redundancy property of the interaction kernel. In fact, the 
q dependence in E^ E (q) and EQ E (q) is very weak for the wide range of q and the 
difference between EQ E {q) and E\ E (q) is also very small. 

Some simplification takes place in the RGM calculation using partial-wave com- 
ponents of the invariant interaction. We can easily show that for each isospin com- 
ponent the direct term (0D+ type) and the knock-on term (1D_ type) give the same 
contribution, when the central (17 = 0) and spin-spin (i? = ss) components are 
added up. For the central IE-, IS- and 15'-type spatial integrals, only the S*-wave 
component of the two-nucleon interaction contributes to the exchange kernel. Thus, 
by incorporating the spin-isospin factors given in Table HI we find for the central na 
RGM kernel (we remove the tilde of M\ E using the approximation in Eq. (|2-24|) ) 

V KGM (e) = (Vb + G Y ) + G K + eK , 
i 

V D + G v = 2 £(2J + 1) (M % + + M°| - M°| - M°|,) , (2-25) 

7=0 

where fil = 0/ spatial functions are calculated from the g$ (central) component 
in Eq. (|2- 15[) . Similarly, the LS RGM kernel is obtained only from the 0Z? + -type 
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contribution using the /ig component in Eq. (|2T5p . The internal energy of the a 
cluster is given by 

E ' M = 3 • W + 4 \- 3E oe 01 ^) + E oe%)] + 2e 2 ^ , (2-26) 

where E^^ 1 (q) is the central component of EQ E (q) with an isospin 1 = 0, 1. 

To calculate the RGM kernel, it is sometimes convenient to use either of the 
forms g{p,p';K) and u(k',q';K), which are equal to each other with k' = p — p' 
and q' = (l/2)(p + p'). The kernel expressions become sometimes simple, if we 
use k = — q i and q = (1/2) (q^ + gj instead of q^ and q^ For example, the 
exchange normalization kernel M^(qj,q i ) = K(qj,q i ) and the exchange kinetic- 
energy kernel G K (qj, are given in Eqs. (jA- 1[) and (|A-2p . The -E-type integrals are 
given in Eqs. (|2T9p - (j2-21|) . For the direct term, we obtain 

M 0D+ (lMi) = ^ (3^) 2 jdq' exp |-A (q' - \q 

Ku(k,q';2\q-q'\) , (2-27) 
and the S- and S"-type integrals are given by 



3 

I / 3 \ 2 

Af| is |(q / ,qJ = M JV (q / ,g i ) ^3 J 

■/ ■/ 1 - A(p + V 2 ) + y + - & (!<?/ + ft) 



15' / 



5 ( P , P ' 2| g/ -p|) I 

g(p,p ,2\qt -p \) J 

2.3. Partial-wave expansion 

The partial-wave decomposition of the na RGM kernels derived in the preceding 
section can be carried out for two different types of the partial-wave decomposition 
of the invariant G-matrix NN interaction: 

00 

g(p,p';K)=J2W+l)9e(p,p';K)P i (p-p') , 

1=0 

00 

u(k', q'; K) = ^(2A + 1) U\(k', q'; K) P x (k ■ q') , (2-29) 
A=0 

for the central part and 

00 

h(p, p'; K) = ]T(2£ + !) hiip, p'\ K) Pj{p ■ p') , 

e=i 

00 

(k',q';K)=Y J (^ + l)u L x S (k',q';K)P 1 x (k -q') , (2-30) 



u 

A=l 
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for the LS part. Here, the it-type partial-wave components are calculated from 

gL,M;X) + i \ i[25(5 + i)-3l \ g «#<p>p>*) 



gl sl M;K) J 4^V2^+iy { J [25(5 + 1) -3] 

(£ = 0, 1, 2 r 

1 ^yG^^p'; if) + *^G^W; K) 



21 + 3 ^ie+i 
{£ + l){2l + Y\ u ' u 



G[R}{p,p';K) 



{1=1, 2, 3,.--) , (2-31) 



using the transformation formula (see Eq. (2.22) of Ref. 3)) 



u \ 



(k,q;K)=J2 W+l)Fi? j Xe (k,q)9t(Pi,P j ;K) {Q = C, LS) . (2-32) 

i,j ^=0 or 1 



Here, gf[pi,Pj',K) is assigned to gQg(p,p';K) or g* s g(p,p'; K) in Eq. (|2-3ip for the 
central part and to h,Qg(p,p f ; K) for the LS part. In Eq. (|2 31|) . the prime on the 
summation symbol indicates that the sum is taken only for quantum numbers such 
as J and 5 that satisfy the generalized Pauli principle (— iy+ s+I = — 1 ; and the 
angular-momentum coupling (IS) J. Similarly, the LS component involves only odd 
I for I = 1 and even t for I = 0. As discussed in Ref. 3), u\(2q', k/2; K) needed 
in the calculation of Mm_(k,q) = Mi£>_(q^,q i ) is obtained by simply changing 
giiPiP'\ K) to (—l) e ge(p,p'', K) in Eq. (|2 32|) . This is a consequence of the symmetry 
property satisfied by the coefficients, F l Q j M (k,q), given in Eq. (2.23) of Ref. 3). Be- 
cause of this symmetry, we only need to calculate the partial wave components for 
Mo_D+(qj,gj) for the spin-independent central component g^ipip ']K), as shown in 
Eq. (12^231) . 

For the E-type RGM kernel, only the 5-wave components of the G-matrix in- 
teraction contribute to the RGM kernel: 



3 

S(K) = - I — 

7T \ TXV 



3 

^ 2 J p 2 dpp' 2 dp' exp j-i-(p 2 +p' 2 ) j g (p,p';K) , (2-33) 



and 



E v E (q) = 4vr 2 K 2 dKe-MK-h) 2 i (J^qKj £{K) , 

3 

EY E (q) = 4vr (JL) 2 jH ' K 2 dKe-UK-bf 1 (l qK \ 8{K) , (2-34) 

where i\(x) = e~ x i\(x) = e~ x i x j\(—ix) is employed. For the 0L> + -type RGM kernel, 
the partial wave decomposition is carried out, including the momentum dependence 
of K = 2\q-q'\: 

oo 

M 0D+ (k,q) = ^2(2X + 1) M 0D+X (k,q) P X (k-q) ■ (2-35) 
A=0 
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The result is 

3 j 

M 0D+X (k,q)= 4vre-^ fc2 (JLY jT ' q' 2 dq' \ J dx e~^-^)' 

x e-iw'C 1 -^) UA (k,q / ;2y/q2 + q' 2 -2qq'x} P x {x) . (2-36) 

This term (and also M^g term of Eq. (|2-24|) ) is transformed back to the partial wave 
component Mqd + t(p f , Pi) by the inverse transformation of Eq. (|2-32p . The partial- 
wave components of the LS term MQ^ +£ (pf,pi) are similarly obtained. For the S 
type (and also for S' type), it is easiest to calculate Misi(qf,qi) directly from 

oo 

M hS (q f ,q i ) = Y,(^+ 1 ) M 1S i(q f ,qi) P t {q f ■%) . (2-37) 

1=0 

We obtain 

3 

M lse (q f , qi ) = (-l)^(iy e-i^/W) f p «dp>e-^ 

x J™ K 2 dK e ~i; K2+ ^f + ^ K ~i e (^qiA 

xi J^dx e -hu^-) go Uq f 2 + KyA-q f Kx,p';K \ P e (x) . (2-38) 

Here again, we only find the contribution from the S'-wave components of the G- 
matrix interaction. If we modify \/q/ 2 + K 2 /4 — q$Kx — > qf in the third line of 
Eq. (|2-38p . it becomes i£ {qfK/2v) go(qf,p';K). The S"-type component is obtained 
from the symmetry discussion as 

Mis>e(qf,qi) = M 1S e(qi,q f ) . (2-39) 

With all of these contributions, the partial-wave component of the RGM Born 
kernel for the na system is given by 

V e RGMJ (q f ,q i ;e) = V e c (q f ,q i ;e)+V i LS (q f ,q i ) (£-S) je , 

Vf(qf, 9»; e) = Gfilf, Qi) + eKtfaf, 1i) 
l 

+2 £(21 + 1) [M§£ +e {q f , 5i ) + MU e (Qf, ft) ~ <Si(df, Qi) ~ M^ e (q u q f )] , 
7=0 

1 

V f LS (q f ,q t ) = 4^(21+ l)M^l £ (q f , qi ) , (2-40) 
1=0 

where (i-S) je = 1/2 for J = 1+ 1/2 and -{i + l)/2 for J = £-1/2, and K e (q f , qi ) 
and Gf(qj,qi) are given in Eq. (|A-3p . 
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2.4. Selection of starting energies of G-matrix 

In previous calculations of B%a interactions with B$ = A, £ and S, we used 
the starting-energy dependence of the G-matrix according to the rule of the nuclear 
matter calculation. Namely, we used the angle-averaging procedure over the relative 
momentum q (which is the integral variable) in the relationship, K = q l + q 2 
and q = (l/2)(q 1 — q 2 ), under the constraint \q 2 \ < kp. This procedure yields 
K = 2[ qi 2 + q 2 - q iq (l + [— l|a=| 1])] 1/2 and q 2 = v / K 2 /2 + 2q 2 - q x 2 . Here, x is 
defined by x = (q\ 2 + 4g 2 — kp 2 ) / '(4qxq) and [— l|x|l] = max{— 1, mm{x, 1}}. The 
starting-energy uj is then calculated from 

w = ^r-q 2 + U N ( qi ) + U N (q 2 ) , (2-41) 
M N 

after the subtraction of the conserved cm. energy (h 2 /2Mn)K 2 . The assignments q 
to qi and q' to q in Eq. (|2-27p (and similar prescription in the IS- and IS" terms in 
Eq. (|2-28p ) clearly destroy the redundancy property of the na RGM kernel. Although 
the K dependence is explicitly treated in the partial wave decomposition, we need 
at least one more vector to uniquely specify the momenta q\ and q 2 even under 
the assumption of the angle-averaging procedure. Here we use a set of K and q 
to specify the starting energy uj. This method, however, involves some ambiguity, 
since the definitions of q x and q differ for each type. The explicit expressions in 
Eqs. (|2-27p - (|2-28p indicate that the roles of q and q' in the 0D + and ID- types are 
taken over by qj and p in the 15 type and by q^ and p' in the IS" type. Such a choice 
of three different sets of the G-matrix interaction parameters apparently causes a 
problem. Here, we choose the 0D + -type definitions for q l and q as a standard set. 
Namely, we assume that Eq. (|2-18p is explicitly given by 

( Pl ,p 2 \u\p[,p' 2 ) = 5(K - K')j^g {p,p'-,u; {\(p + p')/2\,K) ,K,k F ) . (2-42) 

In practical calculation, we first assume q = \q\ and K = \K\ and calculate the G- 
matrix using the angle- averaging procedure of K under the constraint \K /2 — q\ < 
kp- Then the rule of the partial- wave decompositions in g and u is given by 

g (p, p'; lo (q, K) , K, k F ) = ^(2£ + 1) g e (p,p'; u(q, K),K, k F ) P e (p, p') 

l 

= u(k',q';oj(q,K),K, kp) 

= Y / ( 2 ^ + l )Mk',q';u;(q,K),K,kp)P x (k',q r ) . (2-43) 
A 

We pick up the q = q' portion in the last line of Eq. (|2-43p and define 

u x (k',q';K,k F )= Ux (k',q';u;(q',K),K,kp) . (2-44) 
The full G-matrix interaction is constructed using 

u(k', q'; K, k F ) = ^(2X + 1) u x {k\ q'; K, kp) P x (k , q') 
A 

= g(p,p';K,k F ) =Y,(M + l)9e(p,p'\K,k F )P e (p,p') , (2-45) 

t 
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which is equivalent to 

g(p,p';uj(q',K),K,k F ) = u(k' ,q';u;(q' , K), K,k F ) . (2-46) 

Since the q' dependence in the starting-energy part is absorbed in the relative mo- 
menta in the bra and ket sides of the G-matrix, we can apply the previous whole 
formalism to the G-matrix interactions u and g. Note that we only need the S-wave 
component for g, since only this component is required for the OE, IE, IS and IS" 
types. 

The angle average over K in the G-matrix calculation is carried out similarly to 
that in the previous case. We start from 

u(q, K) = ^-q 2 + U N (K/2 + q) + U N (K/2 - q) , (2-47) 

M N 

and make the angle-averaging over K under the constraint \K/2 — q\ < kp. Then 
we find 

K ) = ^1 2 + U N ^K 2 /A + q^ + Kq(x}) 

+U N (\f K 2 /A + q 2 - Kq(x}) , (2-48) 



where 



( x ) = l(l + [-l\ x \l]) 




with 



r i 

|(1 + x) for { -l<x<\ , (2-49) 



' n K 2 +q 2 -k F 2 ) . (2-50) 



Kq \A 

In two special cases, we find 



^ K) ={ ^Q 2 + U N (K/2 + q) + U N {\K/2 - q\) 

{ K/2 + q<k F 
ioi I . (2-51) 

[ \K/2-q\ > k F 

§3. Results and discussion 

3.1. NN G-matrix and na phase shifts 

Figure Q] shows the nucleon s.p. potential of the model fss2, obtained by G-matrix 
calculation for symmetric nuclear matter. The Fermi momentum kp = 1.20 fm -1 
is assumed in the continuous prescription for intermediate spectra. The solid curve 
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Fig. 1. Comparison of the nucleon s.p. potentials predicted by the renormalized RGM kernel and 
the eK prescription. The model is fss2 and &f = 1.20 fm -1 is assumed in the continuous 
prescription for intermediate spectra. 

(ren. RGM) is obtained from the renormalized RGM kernel and the dashed curve 
(eK) from the energy-dependent RGM kernel with the explicit eK term. We find 
that the off-shell transformation by 1/yN gives a rather minor modification for the 
NN G-matrix, giving a slightly repulsive effect to the nucleon s.p. potential. We use 
the prescription of the renormalized RGM to calculate the na Born kernel in this 
paper. The internal energy of the a cluster, calculated from Eq. (|2-26j) with q = 
and v = 0.257 fm" 2 , is E a (0) = -26.5 MeV. 

Figure [2] shows the na phase shifts of some low partial waves, predicted by fss2, 
for the neutron incident energies, E n = 0-30 MeV. The G-matrix calculation in 
symmetric nuclear matter is carried out using the s.p. potential in Fig. Q] and the 
prescription for the starting energies discussed in Sec. 2.4. The Fermi momentum 
kp = 1.20 fm -1 corresponds to 70% of the normal saturation density (0.7 po). The 
h.o. width parameter v = 0.257 fm~ 2 used for the (0s) 4 a cluster reproduces the 
root-mean-square (rms) radius of the a particle. Phase-shift solutions are obtained 
by solving the Lippmann-Schwinger equation for the na Born kernel in Eq. (|2-40p . In 
Fig. El^a), we find that the na phase shifts in Lj = S1/2, P3/2 an d -P1/2 states are well 
reproduced, although the S\/2 phase shift might be slightly too attractive at higher 
energies. For higher partial waves with D 5 / 2 ,3/2 an d -^7/2,5/2 in Figs. [2)[b) and (c), 
the central attraction is not sufficiently large, although the spin-orbit splitting seems 
to be reasonably good. The -D3/2 phase shift is particularly too low in Fig. [2Kb) (and 
also in Fig. Efb)). Part of the reason for the missing attraction in the -D3/2 channel 
is because the specific distortion effect of the d + 3 H channel is not included in the 
calculation, whose threshold opens at E n ~ —22 MeV. 

For comparison, we show in Figs. E^a) - (c) the na phase shifts obtained by 
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Fig. 2. na RGM phase shifts predicted by 
the quark-model G-matrix interaction by 
fss2 for k F = 1.20 fin -1 . (a) S 1/2 , 
P3/2 and P1/2 states. (b) D 5 / 2 and 
D 3 / 2 states, (c) F7/2 and F 5 / 2 states. 
The (0s) 4 shell-model wave function with 
the h.o. size parameter v = 0.257 fm 
is used for the a-cluster wave function. 
The experimental data are taken from 
Refs. 21) and 22). 
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Fig. 3. na RGM phase shifts predicted by 
the Minnesota 3-range force with u = 
0.94687 and v = 0.257 fm" 2 . The LS 
force is KKNN xl.5. (See the text.) 



the standard RGM calculation, using an effective NN force. In this calculation, 
the Minnesota three-range force (MN3R) 18 ) with the Maj orana exchange mixture 
u = 0.94687 is used for the central force. The act phase shifts are well reproduced 
in the aa RGM calculation using this effective NN force. 19 ) For the LS force, 
the two-range Gaussian LS force by Kanada et al., 20 ) is used with the spin-isospin 
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coefficients W = 0.5 and H = -0.5 (no 3 E LS). If we use the force parameters 
given in Table I of Ref. 20), the LS splitting of the P 3 / 2 and Pu 2 states becomes 
too small. This is because these authors introduced the D-w&ve component of the 
a cluster and an extra contribution to the LS splitting originates from the two- 
nucleon tensor force. If we multiply the strength of this LS force by a factor 1.5, it 
gives the correct magnitude, as seen in Fig. [3][a). We call this set of effective NN 
forces MN3R+1.5KKNN. Although the results in Figs. [2] are not as good as those of 
MN3R+1.5KKNN in Fig. [3l it is clear that our quark-model AW interaction gives 
a reasonable description of the not scattering through the G-matrix approach of the 
bare interaction. 



3.2. Scattering cross sections and polarization 

Since the agreement of the calculated phase shifts with the empirical values deter- 
mined from the i?-matrix analyses 21 -*' 22 ) is not complete, we examine the scattering 
cross sections and polarization of the na scattering directly with the experiment. 
Figure H] shows the na total cross sections up to E n = 30 MeV. The prominent peak 
structure at E n = 1-2 MeV is due to the sharp P 3 / 2 resonance. Although the 
calculated results at around E n = 4-6 MeV are slightly small, the agreement from 
12 to 18 MeV is satisfactory. We also examine the differential cross sections and 
polarization in Figs. [5] and [6] at some available energies. We obtain a fare agreement 
between the calculational and experimental results except in the threshold energy 
region, E n ~ 22 MeV, for the decay to the d + 3 H channel. For higher energies, we 
need to introduce imaginary potentials. 
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Fig. 4. Calculated na total cross sections by fss2, compared with the experiment. The experimental 
data are taken from Ref. 23). 
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Fig. 5. Differential cross sections of 
the na scattering predicted by 
fss2. The experimental data are 
taken from Refs. 24) and 25). 
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Fig. 6. Polarization of the na scat- 
tering predicted by fss2. The ex- 
perimental data are taken from 
Refs. 26) and 27). 
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3.3. Analysis of Born amplitudes 

Since we calculated the na Born kernel using only the real part of the G-matrix, 
the success of the na single-channel RGM calculation is qualitative for energies 
greater than 30 MeV. Here, we analyze the roles of the various exchange terms of 
the RGM central kernels from the G-matrix NN interaction. Such an analysis was 
carried out by Thompson and Tang 7 ) for a simple one-term Gaussian NN interac- 
tion. They found that three different groups of the nuclear exchange terms have 
characteristic behavior of various reaction processes. First, the direct (O.D+-type) 
and knock-on (lD_-type) terms, specified by fr> and f\ terms in the Born amplitudes 
below, respectively, have large amplitudes at the forward angle, 9 = 0, characteristic 
of the Wigner-type potential. Secondly, the heavy-particle pickup (IS and IS") term 
ji and the nucleon-rearrangement amplitude f% are important at backward angles. 
They are mainly responsible for the existence of large scattering cross sections at 
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Fig. 7. na Born amplitudes of the G- 
matrix NN interaction by fss2 with fcf = 
1.20 far 1 and v = 0.257 fm~ 2 at the 
cm. energies E c . m . = 50 and 100 MeV. 
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na Born amplitudes of the MN3R 
force with u = 0.94687 and v = 
0.257 fm~ 2 at the cm. energies E c .m. = 
50 and 100 MeV. 
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these angles. In the G-matrix formalism, these Born amplitudes are given by 

oo 1 

fD + fl = 2 E E( 2£ + l ^ 21 + l ) M 0D + i(q, q) Pe(oos 6) , 

OO 1 

f* = ^2 4 E E( 2£ + X )( 2J + *) M °/^> ?) F ^ cos ) > 

£=0 7=0 

1 

Gf(q, q) + 2 £(2/ + l)Jli?i,(<7, ?) + ^(?, ?) 
7=0 

xP £ (cos 0) . (3-1) 

Here, /U = (4/5)Afjv is the reduced mass, q is the relative momentum determined 
from e = (h? /2fi)q 2 , and 9 is the cm. angle for the na scattering. In Eq. (|3-1|) . the 
partial-wave sum is actually taken up to t = 8. We plot fo + fi, /a, and /3 in Fig. [7] 
for the cm. energies E c _ m _ = 50 and 100 MeV. For the effective local iViV forces, 
these Born amplitudes are analytically calculated as in Eqs. (30) - (34) of Ref. 7). We 
show the results of the MN3R force in Fig. [8] for comparison. In the MN3R force, the 
/d term is a Fourier transform of the momentum-independent local potential and f\ 
is that of the nonlocal RGM kernel with a smaller amplitude. Since the /jj and /i 
terms are not separated for the G-matrix interaction, the sum fo + f\ is compared 
with the predictions by fss2. The amplitudes fi and f% arise as a consequence of the 
antisymmetrization procedure and are closely related to the nucleon-exchange effect 
in the backward angular region. Figures [7] and [8] show that the angular and energy 
dependences of these two different NN interactions are very similar to each other, 
in spite of the strong nonlocality of the G-matrix interaction. 

§4. Summary 

The na system is one of the most successful examples, in which microscopic RGM 
calculations using effective NN forces give a good description of the experimental 
data. It is, therefore, very interesting to examine if this result is still valid when 
more realistic NN interactions based on bare interactions are employed in G-matrix 
formalism. Here, we studied the na system by our previous technique, 3 ) in which 
baryon-octet (Bg) a Born kernels are calculated with explicit treatments of the 
nonlocality and the center-of-mass (cm.) motion between Bq and a. Since the na 
system involves the antisymmetrization of nucleons due to the Pauli principle, we 
needed to extend the previous techniques for the direct and knock-on terms of the 
hyperon a interaction to other nucleon-exchange and interaction types such as the 
S and S' types. We found that the explicit treatment of the Galilean non-invariant 
interaction gives some definite recoil effect to the cm. momentum of two interacting 
nucleons, involved in the G-matrix, for each particular interaction type of the RGM 
kernels. If one uses the invariant G-matrix as the input iV./V interaction for the 
na RGM kernel, both of the direct potential and the knock-on interaction kernel 
become nonlocal, and give the same contributions for each isospin component of the 
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NN interaction with I = or 1. This is a common feature of two-cluster systems 
composed of a single nucleon and a nucleus. 

In principle, we can deal with the momentum dependence and the starting- 
energy dependence of the G-matrix according to the explicit expressions derived in 
this study. The treatment of the Fermi momentum kp and the starting energy lu, 
however, needs careful treatment. Local density approximation, usually assumed in 
heavier systems, cannot be justified for (Os)-shell nuclei like the a particle. The 
single-particle (s.p.) energy used in the definition of the starting energy is not well 
defined in the na cluster model. We assumed a constant fcj? to generate G-matrix 
NN interaction for na scattering. The starting energy ui is determined from the local 
momentum q and the cm. momentum K of the two interacting nucleons using angle- 
averaging procedure over K under the constraint \K/2 — q\ < kp. The G-matrix 
calculation is carried out in the continuous prescription for intermediate spectra 
using the energy-independent version of the quark-model NN interaction fss2. We 
found that, in the present framework, kp = 1.20 fm _1 is a most favorable choice 
for generating an appropriate strength of the na interaction. A larger kp say kp = 
1.35 fin -1 , gives less attractive na interaction, and a smaller kp gives more attractive 
na interaction. We carried out the angular-momentum projection of the G-matrix 
NN interaction explicitly including the K dependence. Owing to this procedure, the 
existence of the Pauli-forbidden (0s) state for the 5-wave relative motion is strictly 
preserved, and the admixture of the redundant component of the RGM formalism is 
completely eliminated. 

With these treatments of G-matrix parameters, we found that the central and 
spin-orbit components of the na Born kernel have reasonable strengths under the 
assumption of a rigid translationally invariant (Os) 4 shell-model wave function of the 
a cluster. The na phase shifts in the energy region, E n < 30 MeV, are reasonably re- 
produced for the 5*1/2; -P3/2 an d P1/2 states, while the central attraction is somewhat 
too weak for higher partial waves. The direct comparison of the differential cross 
sections and polarization with experiment shows that these higher partial waves do 
not markedly impair the fit to the experimental data, except for the energy region 
where the d + 3 H channel opens. In the higher energy region of up to 100 MeV in the 
cm. system, we compared the Born amplitudes from the G-matrix NN interaction 
with those from an effective NN force, the Minnesota three-range force. We found 
that characteristic behaviors of three different groups of exchange terms, the direct 
and knock-on terms, heavy-particle pickup terms, and nucleon-rearrangement terms, 
are essentially the same between these two approaches. 

Note that the appropriate strength of the na central attraction in the present cal- 
culation depends largely on how we deal with the strong starting-energy dependence 
in the G-matrix calculation. The present procedure is one of the possible procedures 
for preserving the redundancy property of RGM formalism and still dealing with the 
recoil effect of the a cluster explicitly in the G-matrix NN interaction. Much simpler 
treatments are indeed possible, by assuming some appropriate relative momentum 
q and the cm. momentum K in the starting energy co(q,K). We can find the most 
appropriate kp, which strongly correlates with the approximate reproduction of a- 
cluster internal energy. In all of such calculations, however, the strength of the na 
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LS potential is rather stable and the necessary spin-orbit splitting between the P3/2 
and P1/2 resonances is always reproduced correctly. 

From the present study, we learn that our new folding procedure of the G-matrix 
elements using simple shell-model wave functions works reasonably well for deriving 
the characteristic features of the na interaction. The purpose here is not to study the 
validity or foundations of the model assumptions of the five-nucleon system. There 
are many unsolved conceptual problems in the realistic description of the five-nucleon 
system based on the fundamental AW interaction. Many recent studies imply the 
importance of the tensor correlations and the effect of three-body forces 15 ) even in 
the problem of spin-orbit splitting discussed in the present paper. For applications of 
the present formalism to hyperon-nucleus potentials, however, such details are still 
beyond the scope of the present study. The basic baryon-baryon interaction itself 
contains a large ambiguity, to which we hope to set some constraints by examining the 
hyperon-nucleus potentials using the present approach. Applications to S'-nucleus 
interactions for light nuclei, for instance, will be discussed in the forthcoming paper. 
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Appendix A 

The na RGM Born kernels for the Gaussian-type effective NN force 



The exchange normalization kernel of the na system in the momentum repre- 
sentation is given by 



=(!)'«*{- 

r 1 ri7 , 2 2 , 



(A-l) 



where v is the h.o. size parameter of the a cluster. In many cases, the transformation 
from qjr and gj to the momentum transfer k and the local momentum q via k = 
qj — q t and q = (gy + gj)/2is convenient for simplifying the expressions for the 
kernels. We also use the notation f(9) for Eq. (|A-ip . although it is a function of 
If = I Qf I) Qi = an d cos = (<?/ ' Qi)- The Born kernel for the exchange kinetic- 
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energy kernel is given by 



3h 2 u 



2 2 1,2 

1 I —q + -k 2 

3v V 3 y 4 



(A-2) 



with Mj\r being the nucleon mass. The partial-wave decomposition of these kernels 
yields 
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(A-3) 



where i\{x) = i x j\(—ix) = e x i\(x). For the Gaussian central NN force in Eq. ([2 
the spin-flavor factors in Eq. (|2-6p are calculated to be 



-Xl.E = X\s = X\si 



Xqd + — ^X d , 
-i(X d + X e ) 



X d = 8W + AB - AH - 2M , X e = 8M + AH - AB - 2W 
By using these factors, Vd and G v are obtained as 
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G v (q f , Qi ) = vo- (X d + X e ) [ f E (9) - f s (9) - f s >(8) ] + v -X e f D _ (9) 



(A-5) 



with 



Me) = f(9) 



x < 



1+M 



exp 



exp 



{l%^(f^ + ^) 2 } 

( -p{i%§<? 2 } 



for T 



(A-6) 
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The S"-type function fs'(&) is obtained from fs(0) by simply replacing k by —k. 
For the Gaussian LS force in Eq. (|2-8|) . the spin-isospin factors are given by 



X$ + = {AW - 2H)S , = {AH - 2W)S , (A- 7) 

where S is the total spin of the na system. The LS Born kernels are given by 

H$*h„ ft) - ^(«- - 2tf) (I) 1 exp {4 (A + I) } JL in . S , 

G 1S ta/.9i) = .<„ i5 (2H'--4ff)/ D _(«) (-^j I;""*. ( A ' 8 > 

where n = [q x fc] = x q^]. For the effective NN force, the direct potentials 
become local in the coordinate representation: 

V^(r) = 4 S (W - 2fT) 2 | exp {"jf^ 2 } . (A-9) 

where the LS potential is V£ (r) £ • S. The internal energy of the a cluster is given 
by 

4 nt = 3 • + v {X d + X e ) (-JL-) 5 + 2e 2 ^ . (A-10) 

Appendix B 

The momentum dependence of the interaction kernels for systems of two 

(Os) -shell clusters 

In this Appendix, we will show the explicit momentum dependence that appears 
in the interaction kernels of the Galilean non-invariant two-nucleon interaction for 
systems of two (Os)-shell clusters. For the two-nucleon interaction in Eq. (|2-18p . the 
spatial part of the interaction kernel defined in Eq. (|2 6[) is given by 

M xT {q f ,q l )=M?{q f ,q l ) 

' f dk'^J-(l + ^]^k' 2 - J-V-ife' 



(2vr) 3 J IV In) Av 2x/7 

— ' U'expj » » aV" 

TTul-a/2nJ J y P | l-a/2^u\ A\i 2^ J 

A 1 1 \3 /• f A 1 1 ,„ „ n2 

ciPexp (P - P ) 2 



A - 2 47rt/ 1-/3/2/// y ^1 A - 2 1 - /3/2/x 4i/ 

x u(fcV;|P|) . (B-l) 
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Here M~? (<Z/,Qj) is the normalization kernel given by 

r.\, - > i 27T 1 V \ 1 A-T 2 , 1+Tl, 3 



with r = 1 — x/ fi , fc = qj — qi and q = (q^-\-q i )/2. We assume the mass numbers 
and ^2 with A\, A2 < 4 for the two clusters, A = yli + ^2 is the total mass number 
and n = A1A2/A is the reduced mass number. Almost all of the coefficients and 
vectors appearing in Eq. (|B-ip are defined in Eq. (A. 14) of Ref. 6). The P integral 
is used to suitably treat the cm. motion of the interacting two nucleons. The new 
parameters P$ and j3 are parametrized as 

I 1 1 A ~ 

P = ^(l_ r 2 )(1 _ a/2/i) P > P= (l_ r 2)(l_ a/2/U )yl-2 /3 ' (B ' 3) 

with Pq and j3 given in Table HH 

It sometimes happens that the formula in Eq. (|B-ip cannot be used because of 
the divergence of the coefficients. For example, all the direct terms with x = and 
many of the A\ = 1 or A2 = 1 cases should be treated separately. For the direct 



Table II. Parameters P and (3 in Eq. p-3 
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terms, we obtain 



M 0Eaa (q f ,q i )=8(k) 

3 

1 1 



1 



and 



1 - 2/A a Attu 
x u(k',q';\P\) , 



M 0D+ (qf,qi) = exp 

3 

1 1 



dP exp 



dk' dq' exp 
1 



I -2/Aa Av 



A, 



(B-4) 



1 



1 



1 



1 - 1/2// 

1 - 1/2// 



2/i / Av 
dP dq' exp 
1 



-k' 



1 - 1/2/i 



1 



(1 - 1/Ai)(l - 1/A 2 ) 4ttz/ 

2 



1 



1 - 1/2/i i/ 



1 

— ( 

2/i 



1 







(<*'-<*)" 


2 1 


(i 


A 2 ) 







(1-1/Ai)(l-1/A 2 )4i/ 
x U (fe,q';|P|) . 

When Ai = 1 or A 2 = 1, Eq. (IB-5D cannot be used either. In this case, we obtain 

h- \ I ' 

2/x J Av 



(B-5) 



M D + (q f ,qi) =exp< - 1 



. 3 
1 \ 2 



1 - 1/2/i 7TZ/ 



for Ai = 1 or A 2 = 1 



(B-6) 



The knock-on term kernel MiD_(qt,q i ) can be obtained from Eqs. ()B-5|) and ()B-6p 
by simply changing u(k',q'; \P\) to u(2q',k'/2; \P\). Finally, 1 — (3/2fi becomes zero 
for the S- and S"-type interaction kernels in the na case. In such a case, we can use 
the limit formula 

3 

i&(;5Fj e " v/ " -'<»>■ (B ' 7> 

The last two lines of Eq. (|B-1|) become simply u(k', q', \Pq\) using this procedure. 
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